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Utilizing the Baym-Kadanoff formalism with the polarization function calculated in the random 
phase approximation, the dynamics of the 1^ = quantum Hall state in bilayer graphene is analyzed. 
Two phases with nonzero energy gap, the ferromagnetic and layer asymmetric ones, are found. The 
phase diagram in the plane (Ao, B), where Ao is a top-bottom gates voltage imbalance, is described. 
It is shown that the energy gaps in these phases scale linearly, AE ^ 10-B[T]K, with magnetic field. 
The comparison of these results with recent experiments in bilayer graphene is presented. 
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I. INTRODUCTION 

The properties of bilayer graphene^"— , consisting of two closely coupled graphene layers, have attracted great 
interest. The possibility of inducing and controlling the energy gap by gates voltage makes bilayer graphene one of 
the most active research areas with very promising applications in electronic devices. Recent experiments in bilayer 
graphene^ii showed the generation of energy gaps in a magnetic field with complete lifting of the eightfold degeneracy 
in the zero energy Landau level, which leads to new quantum Hall states with filling factors = 0, ±1, ±2, ±3. Besides 
that, in suspended bilayer graphene, Ref. reports the observation of an extremely large magnetoresistance in the 
I' = state due to the energy gap AE, which scales linearly with a magnetic field B, AE ~ 3.5 — 10.5i?[T]K, for 
B < lOT. This linear scaling is hard to explain by the standard mechanisms^ of gap generation used in a monolayer 
graphene, which lead to large gaps of the order of the Coulomb energy /I ^ B^^^, I = {he/ eB)^/'^ is the magnetic 
length. 

The theory of the quantum Hall effect (QHE) in bilayer graphene has been studied in Rcfs. [lol - [l^ . In particular, 
the gap equation for the quasiparticle propagator including the polarization screening effects has been recently studied 
m Refs. mil While a polarization function with no magnetic field was used in Ref. [l^ the polarization function 
with a magnetic field was utilized in Ref. [H. 

In this paper, wc study the dynamics of clean bilayer graphene in a magnetic field, with the emphasis on the v ^ Q 
state in the quantum Hall effect (QHE) [a brief description of a part of the results of this study was presented in Ref. 

It will be shown that, as in the case of monolayer graphenei^, the dynamics in the QHE in bilayer graphene is 
described by the coexisting quantum Hall fcrromagnctism^ (QHF) and magnetic catalysis^ (MC) order parameters. 
The essence of the dynamics is an effective reduction by two units of the spatial dimension in the electron-hole pairing 
in the lowest Landau level (LLL) with energy E = 0—"—. As we discuss below, there is however an essential difference 
between the QHE's in these two systems. While the pairing forces in monolayer graphene lead to a relativistic-like 
scaling AE ~ VTeSj for the dynamical gap, in bilayer graphene, such a scaling should take place only for strong 
magnetic fields, B > Bttr, where the threshold magnetic field is estimated as Bthr ^ 30T (see Sec. IIII B|) . For 
B < Bthr, a nonrelativistic-like scaling AE ^ \eB\ is realized in the bilayer. The origin of this phenomenon is very 
different forms of the polarization function in monolayer graphene and bilayer one that in turn is determined by the 
different dispersion relations for quasiparticles in these two systems. 

The polarization function is one of the major players in the QHE in bilayer graphene. As will be shown below, its 
role is important because it is proportional to the large mass of quasiparticles, m ^ lO^^JTig ~ lO^K/c^ ^ fi^/e^l 
unless B > SOT, which leads to strong screening. 

Using the random phase approximation in the analysis of the gap equation, we found two competing solutions: 
I) a ferromagnetic (spin splitting) solution, and II) a layer asymmetric solution, actively discussed in the literature. 
Studying how the energy gaps of these solutions depend on the longitudinal component By of the magnetic field at a 
fixed value of the transverse component wc found that while the gap of the solution I increases with B\\ , the gap 
of the solution II decreases as -By increases. Comparing this behavior with that observed in experiment in Ref. and 
calculating the energy density of the ground states for these solutions, we come to the following scenario. While at 
low magnetic fields, the layer asymmetric solution II is realized with the energy gap AE ^ 10i?[T]K, there exists a 
first order phase transition to the ferromagnetic phase corresponding to the solution I at some critical value Bcr- The 
experiment!^ implies that the value of Bcr satisfies Bcr ^ lOT for B\\ = 0. The phase diagram in the plane {Ao,B), 
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where Aq is a top-bottom gates voltage imbalance, is described. 

The paper is organized as follows. In Sec. [Ill the Hamiltonian of the model, its symmetries, and order parameters 
are described. In Sec. Illli by using the Baym-Kadanoff formalisroi^, the gap equation for the quasiparticle propagator 
including the polarization function is derived and the properties of the polarization function are described. In Sec. 
IIV Al the properties of the solutions of the gap equations and the phase diagram of the model are discussed. In Sec. 
IIVB( we compare our results with experiment. In Sec. |Vl we summarize the main results of the paper. In appendix 
\X[ a detailed derivation of the polarization function in a magnetic field in bilayer graphcne is presented. 

II. MODEL 
A. Hamiltonian 

The free part of the effective low energy Hamiltonian of bilayer graphene isi: 

^0 = 1 d'x^+,{x) (^q)' ) ^vM, (1) 

where tt ~ p^^ + ipx^ and the canonical momentum p — —ifiV + eA/c includes the vector potential A corresponding 
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to the external magnetic field B. Without magnetic field, this Hamiltonian generates the spectrum = 
m = 7i/2u|,, where the Fermi velocity vp ~ c/300 and 71 ~ 0.34 — 0.40eV. The two component spinor field vE'v^ 
carries the valley {V = K, K') and spin (s = +, — ) indices. We will use the standard convention: 4']^^, = {4'Ai,'ipB2)Ks 
whereas 5*^/^ = {iJ^B2,''Pai)k's- Here Ai and B2 correspond to those sublattices in the layers 1 and 2, respectively, 
which, according to Bernal {A2 — Bi) stacking, are relevant for the low energy dynamics. The effective Hamiltonian 
([1]) is valid for magnetic fields IT < B < Bthr- For B < IT, the trigonal warping should be taken into accounlfi. For 
B > Bthr, a monolayer like Hamiltonian with linear dispersion should be used. 

The Zeeman and Coulomb interactions plus a top-bottom gates voltage imbalance Aq in bilayer graphene are 
described as (henceforth we will omit indices V and s in the field ^'y^): 



H,„, = ^^BB Jd^x^+ix)a'^ix) + ^ Jd^xd^x' ^'^^^^^^'^K Ao J d^x^+ ix)^T,<fix) 



(2) 



where hb is the Bohr magneton, a is a spin matrix, k is the dielectric constant, and n(x) = 5{z—^)pi{x)+6{z+^)p2{x) 
is the three dimensional charge density {d ~ 0.35nm is the distance between the two layers). The Pauli matrix in 
the voltage imbalance term acts on layer components, and ^ = ±1 for the valleys K and K' , respectively. 
Integrating over z and z' in this equation, one can rewrite Hint as 

H^nt = pbB fd^x^+{x)a^^{x) + ^ [d^xd^x'[V{x-x'){pi{x)pi{x')+p2ix)p2ix'))+2Vi2{x-x')pi{x)p2{x')] 
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+ Ao y d^x^+{x)^T3^ix). (3) 

Here the potential V{x) describes the intralayer interactions and, therefore, coincides with the bare potential in 
monolayer graphene whose Fourier transform is given by V{k) = 2^6^ /nk. The potential V12 describes the interlayer 
electron interactions. Its Fourier transform is Vi2{k) = (27re^/K)(e^'^''/fc). The two-dimensional charge densities pi{x) 
and P2ix) are: 

Pl{x) = ^+{x)Pi^{x) , P2{x) = ^+{x)P2^{x) , (4) 

where Pi = ^^^'^ and P2 — ^ are projectors on states in the layers 1 and 2, respectively. When the polarization 
effects are taken into account, the potentials V{x) and Vi2{x) are replaced by effective interactions Ve//(x) and 
V12 eff{x), respectively, whose Fourier transforms are given in Eas. (|A3p and (|A4[) in appendix 1X1 

B. Symmetries and order parameters 



The Hamiltonian H = Hq + Hint, with Hq and Hint in Eqs. ([TJ and ([3]), describes the dynamics at the neutral 
point (with no doping). Because of the projectors Pi and P2 in charge densities Q, the symmetry of the Hamiltonian 
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H is essentially lower than the symmetry in monolayer graphene. If both the Zeeman and Aq terms are ignored, 
it is t/(^'(2)s X t/(^')(2)s X z!^p x Z^2y\ where U^^'>(2)s defines the U(2) spin transformations in a fixed valley 
V = K, K\ and describes the valley transformation ^ — ^ for a fixed spin s = ± (recall that in monolayer 
graphene the symmetry would be C/(4)i^). The Zeeman interaction lowers this symmetry down to G2 = x 
U^^^l)- X U^^'^l)+ X U'-^'^l)- X Z^p X Z^~\ where [/(^^(l)^ is the C/(l) transformation for fixed values of both 
valley and spin. Recall that the corresponding symmetry in monolayer graphene is Gi = [/*■"'"■' (2)y x J7^~^(2)y, where 
C/*^'*^(2)y is the U{2) valley transformations for a fixed spin. Including the Aq term lowers the G2 symmetry further 
down to the G2 = ^^^^1)+ x U^^^\l)- x U^'^">{1)+ x 

Although the Gi and G2 symmetries are quite different, it is noticeable that their spontaneous breakdowns can 
be described by the same QHF and MC order parameters. The point is that these Gi and G2 define the same four 
conserved commuting currents whose charge densities (and four corresponding chemical potentials) span the QHF 
order parameters (we use the notations of Ref . [isl ) : 

Us : (^l^".,) = {i^KA^.^KA^s + ^^k'A^s^K'A^s + '^KB^s'^KB^s + ^^K' B^s'^K- B^s) , (5) 
f^s ■■ = {^KA^s^KA^s ~ i^K'A^s^K'A^s + ''pRB^s'^KB^s " ^K' B^s'^K' B^s) ■ (6) 

The order parameter ([5]) is the charge density for a fixed spin whereas the order parameter (jG)) determines the charge- 
density imbalance between the two valleys. The corresponding chemical potentials are and /is, respectively. While 

(s) 

the former order parameter preserves the G2 symmetry, the latter completely breaks its discrete subgroup Z^y. Their 
MC cousins are 

As : (*lr3*s> = i.'>P^KA^s'^KA^s - ^K'A^s^K'A^s - i'KB2S^KB2S + i^^K' B^s^ K' B^s) , (7) 
As: (^J^Tg^'s) = {'<l^^KA^s'^KA^s+^^K'A^si^K'A^s-i^'^KB2S^^KB2S--i^^K'B2S^K'B2s) ■ (8) 

These order parameters can be rewritten in the form of Dirac mass terma^^. The corresponding masses are As and 
As, respectively. While the order parameter ([7]) preserves the G2, it is odd under time reversal^S. On the other hand, 
the order parameter (|5]), connected with the conventional Dirac mass A, determines the charge-density imbalance 
between the two layers^. Like the QHF order parameter this mass term completely breaks the z'^y symmetry and 
is even under T. Let us emphasize that unlike a spontaneous breakdown of continuous symmetries, a spontaneous 
breakdown of the discrete valley symmetry z'^y, with the order parameters (vl'j^^'s) and (^'J^Ta^E's), is not forbidden 
by the Mermin- Wagner theorem at finite temperatures in a planar system^!. 

Note that because of the Zeeman interaction, the SU^^\2)s is explicitly broken, leading to a spin gap. This gap 
could be dynamically strongly enhanced^^. In that case, a quasispontaneous breakdown of the SU^^^ (2)5 takes place. 
The corresponding ferromagnetic phase is described by the chemical potential /i3 ~ {fl-^. — /x_)/2, corresponding to 
the QHF order parameter (^I/Vs^I/), and by the mass A3 = (A+ — A_)/2 corresponding to the MC order parameter 

III. GAP EQUATION 
A. General remarks 

In this section, in the framework of the Baym-Kadanoff formalisroi^, and using the polarization function calculated in 
the random phase approximation (RPA) , we derive and analyze the gap equation for the LLL quasiparticle propagator 
with the order parameters introduced above. Recall that in bilayer graphene, the LLL includes both the n = and 
n ~ 1 Landau levels (LLs), if the Coulomb interaction is ignoredi. Therefore there are sixteen parameters Hsin), 
As(n), fls{n), and As(n) with n = 0, 1. 

As will be shown below, including the polarization function in the description of the LLL dynamics is necessary. 
The point is that this function is proportional to a large mass of quasiparticles, m ~ 10~^me ^ 10*K/c^ ^ fi^/e^l 
unless B > SOT, which leads to strong screening effects. 

It will be shown below that the region of relevant values of wave vector k in the gap equation for the LLL states 
is < J/ = k^P/2 < 1. While at small y ^ 1 the dominant contribution (around 80%) in the polarization function 
comes from the transitions between the LLL and the first higher LL with n=2, the number of the LLs providing 
relevant contributions in this function grows with increasing y (for details, see the analysis in appendix E)) . 

Last but not least, a characteristic scale in the bilayer dynamics in a magnetic field is the cyclotron energy fkUc — 
25.5i?[T]K. The applicability of the LLL approximation for a quasiparticle propagator in the gap equation implies 
that the LLL energy gaps should be smaller than hiOc- As we will sec, this condition is fulfilled in bilayer graphene. 
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B. The analysis of the gap equation 



The effective action in the Baym-Kadanoff forniahsm in two-loop approximation is a functional of the full Green's 
function G and has the form 

r(G) = Tr [LnG^i + S'^iG - 1] d^ud^u i^tv[G{u,u')G{u\u)]Vef f{u ~ u') + iv[PiG{u,u') P2G{u\u)] 

X y7L(^i-liO-^tr[G(^i,^i)]tr[G(M^wO]K//(^^-M')-tr[PlG(^i,^i)]tr[P2GK,^i')]^^/L(^i-^i')| , (9) 

where u= (i,r), t is the time coordinate and r = (x^y), Vil{u) ~ Vi2e//(u) — Veff(u) is the interlayer interaction, 
and the Fourier transforms of Veff{u) and V12 eff{u) are given in Eas. (jA3|) and (jA4p in appendix [K\ Note that while 
here the trace Tr, the logarithm, and the product S^^G are taken in the functional sense, the trace tr runs over spinor 
and spin indices. 

The stationary condition dT{G)/SG = leads to the gap (Schwinger-Dyson) equation in mean field approximation, 
which will be written in the form convenient in the presence of a magnetic field: 



G(mi,U2) = S{ui,U2) + i J d^u'id^u'2S{ui,u[)G{u[,u'2)G{u'2,U2) Veff{u'i — u'2) 

+ i [ d^u'id^U2 S{ui,u[) [Pi G(u'i, U2) P2 + ^^2 G(ui, U2) -Pi ] G(m2, U2) V/l(ui - M2) 

d^u^ Siu,,u'^) { tr [ G(ui, ui) ] VeffiO) + (Pi tr [P2 G(mi, «i) ] 
+ P2tr[PiGK,«i)]) V/l(0)}gK,ii2), (10) 

where Veff{0) and V/l(0) are the Fourier transforms of Veff(u) and Vil{u) taken at w = k = 0. 

We will use the Landau gauge for a two dimensional vector potential, A|| = (0, B±x), where B± is the component 
of the magnetic field B orthogonal to the xy plane of graphenc. Then, the free Green's function S'(ui,W2) can be 
written as a product of a translation invariant part S{ui — U2) times the Schwinger phase factor— i^^, 

. {xi +X2){yi -y2) 
21 



5(mi,U2) = exp ( -I — ]S{ui-U2). (11) 



After extracting the Schwinger phase factor in the full propagator, 

.{xi +X2){yi - j/2) 



G(ui, U2) = exp l^-i — jG(ui-M2), (12) 

and making the Fourier transform with respect to we get the following equation for the translation invariant part 
G: 

G{Vt,v) = 5(fi,r)+z j ^ y"d2r'id2r^e*[(^-4)y'i-fe-?/2)<]/2i'5'(f|^i.-r'i) 

X [G(L^,r'i - r^) VeffiSl - c^,r'i - r^) + (PiG(c^,r'i - r^)P2 + P2G(c^,ri - r^)Pi) ViL{n - io,r'^ - r^) 
X G(17,r'2)-i j d^r'2e'^^y'^-y^'^^/^^'sin,r-r'2)[tT[G{0)]VeffiO) 

>itr[P2G(0)] +P2tr[PiG(0)]) V/L(0)}G(l],r^), (13) 



where r = ri — r2. 

The translation invariant part of the free propagator can be expanded over the LLs (compare with Rcfs. miil): 



„2 ^ °° 



%(r;w) = ^7-72 "^^P ("772 ) 7 — r~J ~7^"^2 \ + f^s + ^tsAq) 

2Trr \ 4r / (w + iOsgncj + ^sj^ — I 



2ml^ ""^ \2P J V i^ + iy) 
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where P± = (1 ± T3)/2, En = y h?uj'^n{n — 1) + Aq, uJc = eB/mc is the cyclotron frequency, and the bare electron 
chemical potential /2s = /iQ — sZ includes the Zeeman energy Z ~ /is-B = 0.67i?[T]K (the conventional chemi- 
cal potential /io, responsible for doping, is included for generality). The functions L^ix) are generalized Laguerre 
polynomials, and by definition, L„(x) = L^ix), ^"2(2^) = L'^i{x) ^ 0. 
For the LLL with n = 0, 1, expression (fT4|) takes a simple form: 



%(r;a;) 



2ttI- 



■ cxp 



2P 



2/2 



where 



ai + /is + ^ Aq + iJsgnw 



(15) 



(16) 



Motivated by expression ()15p for the free propagator in the LLL approximation, we will use the following ansatz 
for the full propagator with the parameters /Xs(n), jlsin), As(n), and As(n) related to the order parameters in Eqs. 



Qs(r;cj) 



1 



■ exp 



2P 



where 



and 



1 



id — E^ns + ii^sgnoj ' 



E, 



-i^lsin) + Asin)) + ^iflsin) - As(n)), n = 0, 1, 



(17) 



(18) 



(19) 



are the energies of the LLL states depending on the order parameters Hs{n),fis{n), As(n), As(n). 

Inserting the ansatz (|17p into Eq. (|13|) and using the orthogonality of the Laguerre polynomials, we derive the 
following system of equations for the functions G'jsn(w): 



'is 
1 

4^ 



du cPk 



A 



(27r)3 



G-\in) 



1 n+e 



2 2 



47r/2 



^1 



(27r)3 
1-e 



[G5so(w) + G4,si(w)k2/V2] e-'^'^'/^V;// (^i - c^, |k|) 

A2^Vil{0), (20) 
[Ge,o(c^)k2/V2 + Gesi(c.)(l - e-'^^'^/^Veff {n ~ |k|) 

^2) Vjl{0). (21) 



Here Ai = s sgn(£'-ris) , A2 ~ J2n s sgn(£'+„s). The second and third terms on right hand sides of Eas. ipO)) 
and ((2T|) describe the Fock and Hartree interactions, respectively. Note that because for the LLL states only the 
component "0^25 (V'^is) of the wave function at the K{K') valley is nonzero, their energies depend only on the eight 
independent combinations of the QHF and MC parameters shown in Eq. ([T5)) . 

As is shown in appendix VK[ neglecting the dependence on d in the function Veff{to.,k) describing the exchange 
interactions, one gets 



2^e' 



(22) 



with n(w,fc2) = nii(aj,k) + ni2(cj,k), where the polarization function n,;j describes electron densities correlations 
on the layers i and j in a magnetic field (see Eqs. (jA5[) . (jASp . and (|A9|) ). As to the Hartree interactions, it is (see 
Eq.dMl): 



1//L(w = 0,fc = 0) 



2'Ke^d 



Keff = K + 27re2d(nn(0) - Hi2(0)). 



(23) 
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It is estimated in appendix \X\ that the value of the dynamical part of Keff, i.e., K^ff — k, varies in the interval 1-^4. 

We utilize the frequency independent order parameters fi, /i, A, A and take the external frequency = in Eqs. (|20p 
and (|2ip . The static approximation for the polarization function will be used, n(a;, k^) 11(0, k^). As a justification 
of the latter, we present the following argument. Let us assume that the main contribution in the integrals over uj in 
Eas- ipO]) and (PT|) comes from the pole terms in the functions G^sn(w) (see Eq. pS)) ). The contribution of each pole 
in the polarization function has the form Il{uj = £'^„s,k^), and the dependence on E^ns enters through (E^ns/huJc)'^, 
where the cyclotron energy huc — 25.5i3[T]K. As will be shown below, the ratio (E^ns/'^c)'^ is small, ^ 0.15, and, 
therefore, can be neglected in Eq. (|A24p . which leads to a static polarization function 11(0, k^). 

It is convenient to rewrite the static polarization 11(0, k^) in the form 11 = (7n/fi,^)n(y), where both II and y = k^/^/2 
are dimensionless. The function Ii{y) was expressed in terms of the sum over all the Landau levels (see Eq. (|A25[) 
in appendix |X| and was analyzed both analytically and numerically. At y <g; 1, II(y) — O.SSy and its derivative 11' 
changes from 0.55 at j/ = to 0.12 at y = 1. At large y it approaches a zero magnetic field value, II(y) — ln4/7r (see 
Fig. [Ipi. 




2 4 6 8 10 12 
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FIG. 1: The static polarization function 'inll(tj). 

Because of the Gaussian factors ' = in Eqs. (|20| and ([2T|) . the relevant region in the integrals in these 
equations is < y < 1. The crucial point in the analysis is that the region where the bare Coulomb term k in the 
denominator of Veff{k) = Veff{Q,k) (|22p dominates is very small, < y < 10~^i3[T].The main reason of that is a 
large mass of quasiparticlcs, m ~ 10~^TOe ^ lO^K/c^ ^ fi^/e^l. The last inequality takes place unless B > SOT. As a 
result, the polarization function term dominates in Ve//(fc) that leads to Ve//(fc) of the form Ve//(fc) = C{y)h^ /mPk'^ . 
The factor h'^/mPk'^ has the same k dependence as the Coulomb potential in two dimensions, and the factor C(y) 
describes its smooth modulations at < y < 1 (see Fig. [T|). It is unlike the case of monolayer graphene where the 
effective interaction is proportional to 1/fc. 

By using the change of variables k — >■ Zk in Eqs. ((20)l and (j21|) . one can see that \eB\ occurs as an overall factor in 
the front of the integrals in these equations. The latter leads to the scaling AE ^ \eB\ for the dynamical energy gap, 
and not AE ^ \/\eB\ taking place in monolayer graphene^i^ii^ (see Sec. IIVAI below). 

As shown in appendix [Aj the contribution of the LLL with n = 0, 1 in the polarization function is identically zero. 
At y ^ 1, the main contribution (around 80%) comes from the transitions between the LLL and the first higher LL 
with n = 2. With increasing y, the number of higher LLs providing relevant contributions in the polarization function 
grows. 

As to the c ondition of the applicability of this low energy model, according to Ref. [H, it is determined by the 
relation hajc\/n{n — 1) < 71 /4. Its left-hand side is nonzero for n > 2 and increases linearly with B. Taking n ~ 2 
and the sign of equality in this relation, we find the threshold magnetic field Bthr = T w 32 T that determines the 
upper limit for the values of B for which the low energy model is applicable. 

With the static polarization function, the integration over the frequency cj in Eqs. ipDl) and (|2ip can be performed 
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explicitly, and we get a system of algebraic equations for the energies E^ns in Eq. ([19)) : 

1 



ios = fiQ- sZ + ^Aq - 
1 



[sgn {Ec„,) h{x) + sgn {E^is) h{x)] 



{A^+A2)Vefm 
1 , 



47r/2 

-E^is = /io - sZ + ^Ao - 
1 



2 2 

sgn (i^^os) /2(a;) + sgn {E^^is) h{x)\ 



iA,+A2)VeffiO) 



LAA, + l±iA,] VjLiO) 



where the quantities Ii{x) are 



dyfiiy)e y 

K^Jxy + A-Kil{y) 



(24) 



(25) 



(26) 



with fi{y) = (1, y, (1 — yY) for i = 1,2,3, respectively. Here the dimensionless variable x = 2h'^/e^m?P ~ 
{Ahwc/ oi^^i){vF / c)"^ — 0.003i?[T], where a ~ 1/137 is the fine-structure constant and we used the values 71 = 0.39eV, 
huj^ = h^/mf = 2.19B[T]meV, vp = 8.0 x lO^m/s (see Ref. [H). 



IV. SOLUTIONS AND PHASE DIAGRAM: THEORY AND EXPERIMENT 



A. Properties of solutions 



In this section the solutions of Eqs. ([241) and \2h\ and the phase diagram of the system these solutions lead to 
will be described. If the Zeeman term is ignored, the equations for parameters with different spin indices coincide. If 
the voltage imbalance Aq term is absent, these equations are also invariant with respect to the permutation of layer 
indices (^ — > — ^) [note that Ai <-> A2 under the change ^ — Clearly, these symmetries of the gap equations 
are due to the SU^^\2)s x SU^^'\2)s x x z'^ symmetry of the bilayer Hamiltonian discussed in Sec. IIIBI 
(note that if the interlayer Coulomb interaction term Vil were absent, we would have the C/(4) symmetry group, as 
in monolayer graphene). 

Due to the Zeeman and Aq terms, these equations are inhomogeneous. It is natural to expect that the lowest 
energy solution will have the sign correlating with the sign of inhomogeneous terms (solutions with different signs are 
degenerate in the case of homogeneous equations). Without loss of generality, we can assume that Aq is positive. 

At the neutrality point (/xg — and Ai -\- A2 = 0), we found two competing solutions of these equations: I) a 
ferromagnetic (spin splitting) solution, and II) a layer asymmetric solution, actively discussed in the literature. The 
energy ()19p of the LLL states of the solution I equals: 

where Fq{x) ^ Ii{x) + l2{x) and Fi{x) = I2{x) + l3{x) with li in Eq. ([26|). The solution exists for Aq < Z+^^Fi{x). 
Since Ai = A2 = in this solution, the Hartree interaction does not contribute in E^^J^. Note that the dynamical 
term {'h? l2mP)Fn{x) in Eq. ([27[) can be rewritten as {h\eB\/2rac)Fn{x), where Fn{x) depends on B logarithmically 
for a; ^ 1. 

The energy ([TO]) of the LLL states of the solution II is different: 



The last term in the parenthesis is the Hartree one, and the solution exists for Aq > ,2 + Z — Tr-nFiix). For 

illustrative purpose, in suspended bilayer graphene, with k ^ 1, we will use Heff = 4 (see Eq. ([23])). 
The energy density of the ground state for these solutions is (a = /, //): 



^ E E E + i-s0.67B + CAo)sgnE^^l] . (29) 
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It is easy to check that for balanced bilayer (Aq = 0) the solution I is favorite. There are two reasons of that: the 
presence of the Zceman term and the capacitor like Hartrec contribution in the energy E'^^^JJ in the solution II. In 
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FIG. 2: The phase diagram in the {Ao,B) plane at _B|| — 0. Here the effective 



= 4. 



Fig. [21 the phase diagram on the plane (Aq, B), at = 0, is presented. The I (II) area is that where the solution 
I (solution II) is favorite. The two dashed lines compose the boundary of the region where the two solutions coexist 
(the solution I does not exist to the right of the (red) dashed line in the region II, while the solution II does not exist 
to the left of the (blue) dashed line in the region I). The black bold line is a line of a phase transition between the 
phases I and II. Because the solutions coexist in the region around that line, the phase transition is a first order one. 
The equation for the critical value Bcr has a simple form, Bcr[T] ~ 0.4Ao[K]. 

It is noticeable that for any fixed value of B (Aq), there are sufficiently large values of Aq (B), at which the solution 
I (solution II) docs not exist at all. It is because a voltage imbalance (Zecman term) tends to destroy the solution I 
(solution II). 




FIG. 3: The LLL energies of the solutions I (left panel) and II (right panel) as functions of B with By = 0. Here Ao — and 
Aq = 5K for solution I and solution II, respectively. 

For Aq = 0, the dependence of the LLL energies E^^J^ of the solution I on B, at = 0, is shown on the left 
panel in Fig. [3] (the LLL states with opposite ^ remain degenerate in this solution). The perfectly linear form of this 
dependence is evident. Also, the degeneracy between the states of the n = LL and those of the n = 1 LL is removed. 
The energy gap corresponding to the = plateau is AE^^^ = ^ - 14.4B[T]K. 

On the right panel in Fig. |31 the dependence of the LLL energies of the solution II on B, at ~ 0, is shown for 
Aq = 5K. It is also perfectly linear. Unlike the solution I, the LLL degeneracy is now completely removed. As to the 
energy gap corresponding to the = plateau, it is AE^"^ = {e'lII ~ ^ 5K + 9.35[T]K. 

Fig. |4] illustrates how the energy gaps of the two solutions depend on the longitudinal component of the magnetic 
field = Bcosa at a fixed value of the transverse component B± — Bsina. As one can see, while the gap of the 
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solution I increases with By, the gap of the solution II decreases as By increases. These properties of course reflect 
the opposite roles of the Zeeman term in the dynamics of the solutions I and II. 

Thus the results of the analysis of this subsection imply a possibility of the following two scenarios. When the 
top-bottom gates voltage imbalance Aq = 0, the ferromagnetic phase I is favorite for all values of the magnetic field. 
For nonzero Aq, the phase II is realized for the values magnetic fields up to the critical value Bcr, where a first order 
phase transition to the ferromagnetic phase I takes place. As we will discuss in the next subsection, the experiment 
in Ref. clearly prefers the second scenario. 

What can be the origin of Aq? As was pointed out in Ref. |6i; it could be generated due to disorder-induced differences 
in carrier density between the top and bottom layers. A more interesting possibility is that a relatively small Aq is a 
dynamical parameter corresponding to spontaneous breakdown of the discrete valley symmetry Z2'y x .^21/'' bilayer 
graphene with no magnetic field. 




02468 10 02468 10 



B^[T] B^[T] 

FIG. 4: The dependence of the energy gaps of solutions I (left panel) and II (right panel) on the field _Bj_ for different angles. 
The parameter Ao = 5K for both solutions I,II. 



B. Comparison with experiment 

The first experiments in bilayer graphene in a magnetic field^i^ revealed quantum Hall states with the filling factor 
ly = ±4n, 71 = 1, 2... predicted in the framework of the one electron problem in Ref. [l|. No traces of lifting the eightfold 
degeneracy of the LLL and the fourfold degeneracy of higher LLs were observed. 

Recent experiments in bilayer graphene^ showed the generation of energy gaps in a magnetic field resulting in 
complete lifting the eightfold degeneracy in the LLL, which leads to new quantum Hall states with filling factors 
ly = 0, ±1,±2, ±3. While in Ref. lasuspended bilayer graphene was used, bilayer graphene samples deposited on 
Si02/Si substrates were used in Ref. |7|- Because suspended bilayer graphene is much cleaner than that on a substrate, 
the new quantum Hall states in the former start to develop at essentially smaller magnetic fields than in the latter. 
Also, the energy gaps corresponding to these states are essentially larger in suspended samples than in those on 
substrates. Both these experiments clearly showed that the ly = state is an insulating one. 

Since in this paper the dynamics of the v = state in clean bilayer graphene is analyzed, it would be appropriate 
to compare our results with those in suspended graphene in more detail. The central results concerning the v = 
state in Ref. are: a) the observation of an extremely large magnetoresistance in the ly ~ state due to the energy 
gap AE, which scales linearly with a magnetic field B, AE ^ 3.5 — 10.5i?j^[T]K at least for B± < lOT, and b) at 
fixed B±, an increase in the parallel component of the field reduces the magnetoresistance at least for B± < 6T. This 
can be interpreted as reducing the energy gap AE with increasing By . 

As to result a), the agreement of the expressions for both the gaps AE^ and AE^^ derived in Sec. IIV Al with the 
gap AE observed in Ref. is satisfactory. Concerning the result b) , it suggests that the longitudinal magnetic field 
suppresses the energy gap. This fact excludes the ferromagnetic phase as a candidate for the description of the clean 
bilayer graphene at least for B± < 6T (see left panel in Fig. On the other hand, the solution II, describing the 
layer asymmetric phase, is a viable candidate for this role (see right panel in Fig. 2]). 

This conclusion together with the phase diagram in Fig. [T] suggest the following picture. At B± < Bcr ^ lOT, the 
layer asymmetric phase (solution II) is realized. At B± = Bcr, a phase transition to the ferromagnetic phase (solution 
I) takes place. Because these solutions coexist at B± < Bcr, one should expect that it is a first order phase transition. 
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Taking literally the relation Bcr[T] — 0.4Ao[K] derived from the phase diagram in Fig. [T]in Sec. IIV Al we find that 
Bcr ~ lOT corresponds to Ag ^ 25K. However, because the existence of relevant dynamical contributions beyond the 
random phase approximation is quite possible, one should consider this relation just as a qualitative estimate. 



V. CONCLUSION 



The dynamics of bilayer graphene in a magnetic field B < Bthr is characterized by a very strong screening of 
the Coulomb interaction that relates to the presence of a large mass m in the nonrelativistic-like dispersion relation 
for quasiparticles. The functional dependence of the gap on B derived in Sec. IIV Al agrees with that obtained very 
recently in experiment in Ref. 0- The existence of the first order phase transition between the layer asymmetric phase 
and the ferromagnetic one in the (Aq, B) plane is predicted. 

There are still many open issues in this dynamics. In particular: 

a) It would be important to include the chemical potential /iq in the present analysis in order to describe the higher, 
v = 1,2, and 3, LLL plateaus^. 

b) The present ansatz with the sixteen order parameters is the minimal one for describing the breakdown of the 

U''^\2)s X U^^ •'(2)s X -^21/'' ^ ^2V symmetry in bilayer graphene. It could be extended in order to look for other 
solutions of the gap equation. A natural extension would be to include order parameters that mix the n = and 
n = 1 LLL states. 

c) Although in Sec. IIII Bl we presented arguments showing that the static limit for the polarization function is 
at least reasonable, it would be important to check this conclusion directly by analyzing the gap equation with a 
non-static polarization function. 

d) It would be interesting to describe explicitly the dynamics around the threshold value Bthn when the crossover 
between the regimes with the nonrelativistic-like scaling AE ~ \eB\ and the relativistic-like one AE - ^/\^\ should 
take place. 

We arc planning to consider these issues elsewhere. 
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Appendix A: Polarization operator of bilayer graphene in a magnetic field 

The polarization function 11^^ describes electron densities correlations on the layers i and j: 

,5(a; + a;')<5(k + k')n.y(w,k) = -i < 0|p,(w,k')Pj(w',k')|0 > . (Al) 

There are two independent polarization functions, IIn = 1122 and II12 ~ 1121. Taking into account the polarization 
effects, the bare interactions transform into 

i + v-n~\v^2eff{k) v^ffik) J' "^-[vuik) v{k) l,ni2(fc) nn(fc) j' ^""'^ 

with 

Veff{ijJ,k) = r-^ ^ -, (A3) 

' ^ [k + ^{Ii,,+Ii,2){l + e-^^)\[k + ^{Il,,-Ii,2){l-e-^'^)]' ^ ' 

2^6^ fce-^-'^-222£ni2(l-e-2'='^) 

Vi2eff(^,k) = —^5 ^ r ^ o (A4) 

^^^ ^ [A;-H2z^(nn+ni2)(l + e-''^'^)] [fc + ^(nn-ni2)(l-e-^-'')] ' ^ ' 

where k= |k|, and since IIn and ni2 depend on w, the effective interactions Ve// and Vyiejf depend on it too. 
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Neglecting the dependence on d (i.e., taking d = 0), we obtain: 

27re2 

Veff{uJ,k) = Vi2e//(a;,fc) : 



where n(a;, fc^) = nii(cj, k) + ni2(w, k) is the polarization function in a magnetic field. On the other hand, 



V7l(w, k) = Vi2e//(w, k) - Veff{uj, fc) = - 



27re' 



1 -e 



-kd 



(A5) 



(A6) 



K A:+^(nn(w,fc)-ni2(w,fc))(l-e-'=<i) 
therefore, since the interlayer term yri,(a;, fc) appears in gap equations pop and pip only at w = fc = 0, we find that 

2TTe^d 



V/L(c^ = 0,fc = 0) 



Keff = K + 27re^d(nn(0) - ni2(0)). 



By definition, the polarization functions ITn and 1X12 are defined as 

PiS*(a;',k)Pi^(a; + w',p + k) 



, , , doj'd^k 
nii(a;,p) = I I tr 



(27r)3 



(A7) 



(A8) 



dbj'd^k 
ni2(w,p) =1 / -^^tr 



Pi5(w',k)P2 5(w + w',p + k) 



(A9) 



where Pi = (1 + £,t^)/2 and P2 = (1 — C''3)/2 are projectors on layers 1 and 2, respectively, the trace includes the 
summation both over the valley index ^ and spin, and S{uj, k) is the Fourier transform of the translation invariant 
part of the free fermion propagator (|14p in a magnetic field. 

We are interested in calculating the polarization function n(w, fc^) in the random phase approximation at the neutral 
point (^0 = 0). Its expression in configuration space is 



^tr PiS{lo',y)S{uj+u',-y) 

where a small Zeeman term in the fermion propagator will be ignored. Then 

2i ^ r duj' 1 



(AlO) 



n(a;,p) = J d^re-'P^Uiuj^r 



E 



d^r e-rV2i^-.pr |(^. + Ao)(.. + c.' + Ao) [L„{ry2l^)Lrniry2P) + L„_2(rV2?')L™-2(rV2/')] 



2h^r 



4^4 



:Li_,ir^/2l^)Li,_,ir^/2l^) 



{2ml^Y 

Integrating over the angle and making the change of the variable = 2Pt, we get 

00 



(All) 



n(c^,p) 



E 



1 



dte-*Jo{\/2pHH) 



[{Uj' + Ao)(w + + Ao) [L„{t)Lmit) + Ln-2it)Lra-2it)] + 2ujlf Ll_^{t) Ll^_^{t)^ , (A12) 



where Ji, is a Bessel function and En ~ y h?uj1n{n — 1) + A§ (compare with Eq. (HU). 

In order to evaluate the t integral with the first term in the curl brackets in Eq. ()Al2p . we will use the formula 
7.422.2 in Ref. [H: 



dxx''+^ e-"" J,{bx)L''^%ax')L''n{ax') = (-l)"+"(2a)-"-i6'^ exp ( ) L^r"+" I ^ I Lr-"+™-" 



4a 



4q! 



4a 



(A13) 
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Taking i/ = cr = 0, « = !,& = 2^^, y = p^P /2 in this expression, we obtain 

dte-'U2^t)L^{t)L^{t) = (-l)™+"e-yi:r"(2/)i"""(y) = (-l)™+"e-^/n™(y) (A14) 



with 

InM = L:~"^{y)Ll^-'\y). (A15) 

At small y, we find 

Inmiy) Snm - y [^nSnrn + (m + l)Sn.vi+l + (n + l)5„i^n+l] , V ^ 0- (A16) 

In order to evaluate the t integral with the second term in the curl brackets in Eq. (jA12p . 

'dtt^e-'M2V¥t)Ll{t)Ll{t) ^ {-ir+-e-yClSy), (A17) 
we set = 0,(7 = 2,6 = 2^ in Ea. (|JT3| . 

dxxe-^ M2xW)L;^\x^)LI{x') = e-^L^:'"+"(y)L-2+™-" (y), (A18) 



and use the following identity for the Laguerre polynomials on the left hand side of this equation, 

if (x) = i-xy'^'^^^L-^^^ix), l>0,k + l>0. (A19) 

Then we arrive at 

C30 

J dtfe-'M2V¥)Ll.2mlit) = (-l)™+"^_!^e-^L2+"-"(y)L-2+™-(y), (A20) 



and, therefore, 

li'liv) ^im + l)(m + 2)L-„-;-,iy)Vr-{y). (A21) 

Although the symmetry of Inm{y) under the interchange n m is not explicit, it can be checked by using the identity 
(|A19[) . At small y, we get the following expansion for /^„j(y): 

41iv) - {n+l)in + 2)Snm~y[2in+l)^{n + 2)6r,m 

+ S,,^„,+in{n + l){n + 2) + Srn,n+im{m + l){m + 2)], y ^ 0. (A22) 

Therefore the polarization function (|A12p takes the following form: 



' n,m=0 ■' 



[oj'^-El + iO][{uj + uj'Y-El + iO] 

oo 

X {(w' + Ao)(^ + UJ' + Ao) [Inraiv) + /„-2,m-2(y)] + 2c^c 4-2,m-2 (Z/) } ■ (A23) 

After integrating over w' in this expression, we obtain 



n.m—0 



[Inmiy) + In-2,,n-2iy)] ^ l|^4-2,m-2 (?/) 



(A24) 



with Inmjlnm = for n < Or TO < 0. It is noticeable that the contribution of the LLL (with n,TO = 0,1 and 
Eq = El = Aq) in the polarization function is identically zero. 
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Let us discuss the properties of the static polarization used in the main text in more detail. For the static polarization 
function Il{uj = 0, p) in balanced bilayer graphene, Aq = 0, we get 



n(o,p) 



+ 2 



m 


I' 


oo 
n=2 


(-1)" 


Mn 



E 



(-1) 



^ M„A/,„(A/„ + M™) 

n,m— 2 ^ ^ 



[loniy) ~ hAy)] 



^n(,), 



(A25) 



where M„ = \/n{n — 1). Note that the quasiparticle mass m appears as an overall factor only and does not enter the 
function 11(2/). We checked that the double sum is convergent and in numerical calculation we took the upper limits 
in the sum around Umax, ^max = 250, this is enough to calculate n(j/) up to values y = 12 as it is shown in Fig. [TJ 
Using Eas. (jA16[) . (jA22p . and (jA25p . we find the asymptotics of the static polarization function at y — > 0: 



n(y) 



22/ 



E 



y/riyn + 1 + ^/n - 1) (n + ^/rl?^) 



0.55?/. 



(A26) 



The main contribution in this expression (around 80%) comes from the transitions between the LLL and the first 
higher LL with n = 2 (the term 1/^/2 in brackets). With increasing y, the number of higher LLs providing relevant 
contributions in the polarization function grows. 

In a similar way, one can find the expressions for the two independent polarization functions nii(0, p) and ni2(0, p): 



nii(j/) 
ni2(j/) 



ra e 



E 



(-1) 



n+m 



^ A/„A/„,(M„ + M™) 



{M^Mra (/„-2,r„-2(2/) + /„m (?/))] + 2 ^ 



n=2 



(-1)" 



{hu{y) - h,n{y)\ 



E 



(-1) 



n+rn 



Txt? ^ ^ MnMra (M„ + Mm) 



(2) 



n-2,m-2 



At zero momentum, y = 0, we have 

nn(2/ = 0) = -ni2(2/ = 0) = 



n=2 yjn{n-l) 



(A27) 



(A28) 



As is seen, the quantities nii,ni2 are logarithmically divergent separately, and we introduced cutoff rimax- The 
physical origin of this cutoff is the following. For high energy modes, monolayer like dynamics takes place, whose 
contribution to Ilii{y ~ 0) and Hrziy ~ 0) is strongly suppressed (recall that in monolayer graphene 11(0) = 0, see for 
example Ref. [isl ). Therefore the cutoff Umax can roughly be estimated from the condition of the applicability of the 
low energy effective modeli: hwc^jn{n — 1) < 71 /4, which gives n,nax ~ {-iill2^/2hvFf « 45/B[T] (here the values 
m, 71, and vp are taken from Ref. [l|). 

For the effective dielectric constant Kg// in Eq. (|A7P we thus get: 



r, 2 7 "mux -, 

2me^d \ - 1 
— — > , 

n=2 V"("-l) 



fi2 



0-68 E 



n=2 V'^i"' - 1) 



(A29) 



where the value m « 0.054me was used. With n„iax ~ 45/i?[T], we find that for the range of fields from 0.25T to 
lOT the cutoff Umax varies in the interval 6 -f- 180, and therefore the quantity Kg// — k varies in the interval 1-^4. 
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